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Abstract. We consider a model for interest rates, where the short rate is 
given under the risk-neutral measure by a time-homogenous, one-d imensional 
affine process in the sense of [Duffle, Fili povic. and Schachcrmaverl. We show 
that in such a model yield curves can only be normal, inverse or humped 
(i.e. endowed with a single local maximum). Each case can be characterized 
by simple conditions on the present short rate rt- We give conditions under 
which the short rate process will converge to a limit distribution and describe 
the risk-neutral limit distribution in terms of its cumulant generating function. 
We apply our results to the Vasicek model, the CIR model, a CIR model with 
added jumps and a model of Ornstein-Uhlenbeck type. 



1. Introduction 

We consider a model for the term structure of interest rates, where the short 
rate (rt)t>o is given under the m artingale measure by a one-dimensional cons erva- 
tive affine process in the sense of Duffle. Filipovic. and Schachermaver 2003j . An 



affine short rate process of this type will lead to an exponentially-affine structure 
of zero-coupon bond prices and thus also to an affine term structure of yields and 
forward rates. 

We emphasize here that the definition of Duffie et al. 20031 ] is not limited to diffu- 



sions, but also includes processes with jumps and even with jumps whose intensity 
depends in an affine way on the state of the process itself. The class of models we 
consider naturally includes the Vasicek model, the CIR model and variants of them 



that a re obtained by adding jumps, such as the JCIR-model of Brigo and Mercurio 



20061 Section 22.8]. Since they are the best-known, the two 'classical' models of 



Vasicek and Cox-Ingersoll-Ross will serve as the starting point for our discussion of 
yield curve shapes: 

A common criticism of the (time-homogenous) CIR and the Vasicek model is that 
they are not flexible enough to accommodate more complex shapes of yield curves, 
such as curves with a dip (a local minimum), curves with a dip and a hump, or 
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other shapes that are frequently observed in the markets. Often these shortcomings 
are ex plained by 'too f ew parameters' in the mo del (cf. ICarmona and Tehranchi 
20061 Section 2.3.5] or iBrigo and Mercuriol 20061 Section 3.2]). However if jumps 



are added to the mentioned models, additional parameters (potentially infinitely 
many) are introduced through the jump part, while the model still remains in the 
scope of affine models. It is not clear per se what consequences the introduction of 
jumps will have for the range of attainable yield curves, and this is one question we 
intend to answer in this article. 

Moreover, there seems to be some confusion about what shapes of yield curves 
are actually attainable even in well-studied models like the C IR-model. While 
most sources (including the original paper oflCoxet all Il985ln mention inverse, 
normal and humped shapes, Carmona and Tehranchil [2006L Section 2.3.5] write 
that 'tweaking the par ameters [of the CIR model] can produce yield curves with one 
hump or one dip', and Brigo and Mercuriol 20061 Section 3.2] state that 'some typ- 
ical shapes, like that of an inverted yield curve, may not be reproduced by the [CIR 
or Vasicek] model. ' In our main result, Theorem 13. 9[ we settle this question and 
prove that in any time-homogenous, affine one-factor model the attainable yield 
curves are either inverse, normal or humped. The proof will rely only on tools 
of elementary analysis and on t he characterization of affine processes through the 



generalized Riccati equations of iDuffie et alj [2003 1. 



Another related problem is how the shape of the yield curve is determined by the 
parameters of the model, and also how - when the parameters are fixed - the yield 
curve is determined by the level of the current short rate. We show in Section l4~2l 
that also in this respect the CIR model h as not been com pletely understood and 
discuss a misconception th at originates in Cox et all [1985} and is repeated for ex- 
ample in Rebonatd . 11998 1 . 



In Section 13.31 we provide conditions under which an affine process converges 
to a limit distribution. We also characterize the limit distribution in terms of its 



cumulant generating func tion, extending results of Jurek and Vervaatl 1983| and 
Sato and Yamazatol 1984 ] for OU-type processes to the class of affine processes. 



These results can again be interpreted in the context of interest rates, where they 
can be used to derive the risk-neutral asymptotic distribution of the short rate 
(?*t)i>o as t goes to infinity. 



We conclude our article in Section0]by applying the theoretical results to several 
interest rate models, such as the Vasicek model, the CIR model, the JCIR model 
and an Ornstein-Uhlenbeck-type model. 



2. Preliminaries 

In this section we collect some key results on affine processes from lDuffie et al' 
[2003j . In their article affine processes are defined on the (m + n)-dimensional state 
space R™ x M. n , and we will try to simplify notation where this is possible in the 
one-dimensional ca se. R esults on affine processes with state space R^o can also be 
20011 ] . 

Definition 2.1 (One-dimensional affine process). A time-homogenous Markov pro- 
cess (r t )t>o with state space D = R^ or R an( l its semi-group {Pt)t>o ar e called 



found in Filir 
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U = 



affine, if the characteristic function of its transition kernel pt(x, .), given by 

Pt(x,u)= / e ui p t (x,dg) 

J D 

and defined (at least) on 

'{ueC:Re«<0} if D = R^ , 
{neC:Re« = 0} if D = R , 

is exponentially affine in x. That is, there exist C-valued functions <p(t,u) and 
tp(t,u), defined on R^o x W, such that 

(2.1) Pt(x, u) = exp(<f>(t, u) +xip(t ) u)) for all x £ D, (t, u) £ R^o x W . 

For subsequent results the following regularity condition for (r t )t>o will be 
needed: 

Definition 2.2. An affine process is called regular if it is stochastically continuous 
and the right hand derivatives 

<9 t + 4>{t,u)\ t -Q and ip(t, u)\ t=0 

exist for all u £ U and are continuous at u = 0. 

Definition 2.3. The parameters (a, a, 6, /3, c, 7, m, /i) are called admissible for a 
process with state space R^o if 

a = 0, 

a, 6, c, 7 £ R^ , 
/3 e R, 

m, /i are Levy measures on (0, 00), where m satisfies 
(£ A l)m(d£) < 00 , 

'(0,00) 

and admissible for a process with state space R if 

a, c £ R^o , 
&,/3 61, 
m is a Levy measure on R \ {0} , 
a = 0,7 = 0, /z = . 
Moreover define the truncation functions 

M0 = K ^ and MO 4^ ^ D=R >°, 

\th2 if D = R \0 if D = R 

and finally the functions F(u), R(u) for it £ C as 

(2.2) F(u) = au 2 + bu- c+ [ (e"* - 1 - uh F (£)) m{d£) , 

Jd\{o} 

(2.3) J?(u) = au 2 + (3u - 7 + / (e ui - 1 - uh R (£)) ■ 

Jd\{o} 



Th e next result is a one-dimensional version of the key result of iDuffie et al 
2003| : 



4 



MARTIN KELLER- RESSEL AND THOMAS STEINER 



Theorem 2.4 ( Duffie. Filipovic. and SchachermaveiL Theorem 2.7). Suppose {rt)t>o 
is a one- dimensional regular affine process. Then it is a Feller process. Let A be 
its infinitesimal generator. Then C^°(D) is a core of A, C^(D) C T>(A) and there 
exist some admissible parameters (a,a,b,j3,c,^f,m,fi) such that, for f S C^{D), 

Af{x) = (a + ax)f"{x) + (b + px)f'{x) - (c + jx)f(x)+ 



(/(* + - f{x) - f(x)h F (0) m{aX)+ 

D\{0} 

x [ {f{x + 0- fix) - f{x)h R {0) n(dQ . 
Jd\{o} 



(2.4) 



Moreover (f){t,u) and ip{t,u), defined by (|2.1I) . solve the generalized Riccati equa- 
tions 



(2.5a) 
(2.5b) 



d t <j)(t, u) 
d t ip(t,u) 



F(iP(t,u)), 
R(i>{t,u)) , 



0(0, u) 
0(0, u) 



0, 



Conversely let ia,a,b,/3,c,j,fi,m) be some admissible parameters. Then there ex- 
ists a unique regular affine semigroup {Pt)t>o with infinitesimal generator (|2.4|) . 
and (|2.ip holds with <fi{t,u) and ip{t,u) given by (|2.5[) . 

Closely related to affine processes is the notion of an Ornstein-Uhlenbeck (OU- 
)type process. These processes are of some importance, since they usually offer 
good anal ytic tractab ility and have been studied for longer than affine processes. 
Following Sato! 19991 Chapter 17] an OU-type process (X t )t>o can be defined as 
the solution of the Langevin SDE 

dX t = -\X t dt + dL t , X e M., Xq e M., 

where {L t )t>o is a Levy process, often called background driving Levy process 
(BDLP). In an equivalent definition, an OU-type process is a time-homogenous 
Markov process, whose transition kernel pt(x, .) has the characteristic function 

Pt{x,u)=expi / F{e~ Xs u)ds 



xe~ xt u 



where F (u) is the characteristic exponent of [L t )t>o- From the last equation it is 
immediately seen that every OU-type process is an affine process in the sense of 
Definition 12.11 It is also seen that in the generalized Riccati equations (12. 5ft for an 
OU-type process necessarily R{u) = —\u. Comparing this with (|2.3|) and Defini- 
tion [231 h is seen that any regular affine process with state space K is a process of 
OU-type. The reverse, however is not true, as there also exist OU-type processes 
with state space M^o- We will give an example of such a process in Section l4~4l 



Naturally we will not only be interested in the process (rt)t>o itself, but also in 
its integral J r s ds and in quantities of the type 



(2.6) 



Qt f{x) := E 



exp 



r s ds f{r t ) 



where / is a bounded function on D. The nex t result is an application o f the 
Feynman-Kac formula for Feller semigroups ( cf. Rogers and Williams! 1994L Sec- 
tion III. 19]) and can be found in iDuffie et al.l [200 j ]. It relies on the positivitv of 
{rt)t>o and is therefore only applicable if D 
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Proposition 2.5 (iDuffie. Filipovic. and Schachermaven . Proposition 11.1). Let (rt)t> 



be a one-dimensional, regular affine process with state space R^o- Then the fam- 
ily (Qt)t>o defined by (|2.6p forms a regular, affine semigroup with infinitesimal 
generator 

Bf(x)=Af(x)-xf(x) for all f G C 2 C (D) . 

We will make extensive use of the convexity and continuous differentiability of 
the functions F and R from Definition ^. 31 These properties are established in this 
Lemma: 



Lemma 2.6. If c — 7 = then F, R as defined in Definition \2.3\ have the following 
properties: 

(i) R{0) = and F(0) = 0. 
(ii) R(u) < 00 for all u G (—00, 0]. 

(Hi) If F(u) < 00 on (ci,C2) C R, then F is either strictly convex on (01,02) or 
F(u) = bu for all wel. The same holds for R with b replaced by [3. 

(iv) If F(u) < 00 on (cx,C2) C R, then F is continuously differ entiable on (01,02). 
Also the one-sided derivatives at c\ and 02 are defined but may take the values 
—00 (at ci) and +00 (at c-i). The same holds for R. 

Proof. Property (i) is obvious. If D = M then by Definition 12.31 R(u) — (3u such 
that (ii) follows immediately. If D = IR^o we use the estimate 

(2.7) |e* - 1 - uh R (0\ < \u\ (0(e) A l) , 

for all u G (—oo,0] and £ G R^o, and (ii) follows from ()2.3|1 . For Property (hi) 
note that by the Levy-Khintchinc formula there exists an infinitely divisible random 
variable X, such that F is its cumulant generating function, i.e. F(u) = logE [e uX ~\ 
for u G (ci, 02). Choosing two distinct numbers u, v G (ci, 02), we apply the Cauchy- 
Schwarz inequality to 



F = logE 



e 2 • e 2 



< log JE[e uX ] ■ E[e vX ~* 



F(u) + F{v) 



which shows convexity of F. The inequality is strict unless there exists some c 7^ 
such that e uX = ce vX almost surely. This can only be the case if X is con- 
stant a.s., in which case F is linear. The same argument applies to R. Property 
(iv) follows from the conve xity and from the fact that F and R are analytic on 
{u G C : Re u G (c u c 2 )} (cf. iLukacsl [l96Ct Chapter 7]). □ 



3. Theoretical Results 

We will now use the theory from the last section to calculate bond prices, yields 
and other quantities in an interest rate model where the short rate follows a one- 
dimensional regular affine process (r t ) t >o under the martingale measure. Naturally 
we will also make the assumption that (r«) t>n is conservative , i.e. that pt(x, D) = 1 
for all (t,x) C R^ x D. This implies by iDuffie et"aD [20031 . Proposition 9.1] that 



c = 7 = in Definition 12.31 We will need some additional assumptions which are 
summarized in the following condition: 

Condition 3.1. The one-dimensional affine process (rt)t>o is assumed to be reg- 
ular and conservative. In addition, if the process has state space R, such that by 
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Definition 12.31 R(u) = flu, we require that 

pi/ \ - t c f (1//3, 0] if /3<0, 

(3.1) -F( M ) < 00 i° r au «6 { 

I (—oo,0] else . 

It will be seen that the condition on F is necessary to guarantee e xistence of bond 
prices for all maturities in the term structure model. Bv lSatol 19991 Theorem 25.17] 
we get an equivalent formulation of Condition 13.11 if we replace F(u) < oo by 
Ji£l >;l e"^ m(d£) < oo. Next we define a quantity that will generalize the coefficient 
of mean reversion from OU-type processes: 

Definition 3.2 (quasi- mean- reversion). Given a one-dimensional conservative affinc 
process (r t )t>o, define the quasi-mean-reversion A as the positive solution of 

(3.2) R(-l/\) = l. 

If there is no positive solution we set A = 0. 

Since R is by Lemma |2"1)1 a convex function satisfying R(Q) = 0, it is easy to see 
that (|3.2p can have at most one solution and thus A is well-defined. The name quasi- 
mean-reversion is derived from the fact that if (r t ) t >o is a process of OU-type with 
positive mean-reversion, then R(u) = f3u and the quasi-mean-reversion A = —(3 
is exactly the coefficient of mean reversion of (rt)t>o- When the process (rt)t>o 
satisfies Condition 13. 1[ it is seen that F must be defined at least on (— 1/A, 0]. 
We will encounter several times the condition that A > 0. The next result gives an 
equivalent formulation in terms of (a, /3, fi): 

Proposition 3.3. The quasi-mean reversion A is strictly positive if and only if 
a > °> Id\{o] M£) M(<2£) =°°, or (3- J D \ {0} h R (g) fi(d£) < 0. 

Proof. First note that by Lemma [2.61 R(u) < oo for all u S (— oo,0]. Using the 
estimate (|2.7I) and a dominated convergence argument it is seen from (|2.3p that 

, . R(u) 

(3.3) lim v ' - - 



(3.4) lim R{U) ^ =(3 :=f3- I // /f (i ) n(<H ) 



lD\{0} 

where (3q can also take the value — oo. Suppose now that a > 0. Then by (|3.3p we 
get liniu^-oo R(u) — oo. Since R(0) = and R is continuous it follows that there 
exists a A > such that R(—l/X) = 1. Similarly if a — 0, but /3q < 0, it follows 
from (|3.4p that lim tl _^_ 00 R(u) — oo and thus again that A > 0. 
Conversely, suppose that a = and 0o > 0. Then 

lim R'(u) = lim = /3o > . 

u — > — oo u — > — oo u 

By the convexity of R it follows that R'(u) > for all it G (— oo, 0). Since R(0) = 
this implies that R{u) < for all u S (— oo,0), and consequently that A = 0. □ 

3.1. Bond Prices. We consider now the price P(t, t + x) of a zero-coupon bond 
with time to maturity x, at time t, given by 



P(t,t + x) = E 



exp I — / r s ds 
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The affine structure of {rt)t>o carries over to the bond prices, and we get the 
following result: 

Proposition 3.4. Let the short rate be given by a one- dimensional affine process 
(Xt)t>o satisfying Condition \3.1\ 

Then the bond price P(t, t + x) exists for all t, x > and is given by 

(3.5) P(t,t + x) = exp(A{x) +r t B(x)) 

where A and B solve the generalized Riccati equations 

(3.6a) d x A(x) = F (B(x)) A(0) = , 

(3.6b) d x B(x) = R (B(x)) - 1 B{0) = . 



Proof. If D = R^o the assertion follows directly from Proposition 12.51 by noting 
that P(t,t + x) = Q x 1. 

If D = M then, as discussed after Theoreml2. 4| (rf) f >o is a process of OU-type and 
R(u) has the simple structure R(u) = f3u. Bv TSatol [19991 (17.2) - (17.3)] we obtain 
in this case directly that 



(3.7) E 



t+x 



exp — / r s ds 



exp(Y F(B{s))ds + r t B(x)\ , 



with B{x) = (1 — e l3x )/j3 if /3 ^ and B(x) — —x when f3 = 0. As a function of 
x € K>0: B is continuously decreasing from to 1/(3 if (3 < 0, and from to -co 
if /? > 0. It is therefore seen that the integral on the right side of (|3.7p is finite for 
all x S M^o if and only if F satisfies (|3.1[) , as required by Condition 13.11 □ 



Corollary 3.5. Let (r t ) t >o satisfy Condition \3.1\ and have quasi-mean-reversion 
X. Then the function B(x) from Proposition^^ is strictly decreasing and satisfies 

lim B(x) = -1/A . 

X — >oc 

Proof. The result follows from a qualitative analysis of the autonomous ODE (|3.6b[) . 
Let A > 0. Since R{— 1/A) — 1 = the point x* := —1/A is an critical point of (|3.6bl) . 
By the convexity of R and the fact that R(0) = it follows that R'(x*) < such 
that x» is asymptotically stable, i.e. solutions entering a small enough neighborhood 
of must converge to a;*. Since R(x) — 1 < for x £ (#*, 0] and there is no other 
critical point in (x*,0], we conclude that B(x) - the solution of (|3.6b[) starting at 
- is strictly decreasing and converges to a;*. 

If A = then there is no critical point in (— oo, 0] and R{x) — 1 < for x S (— oo, 0]. 
It follows that B(x) is strictly decreasing and diverges to — oo. □ 

3.2. The Yield Curve and the Forward Rate Curve. The next results are the 
central theoretical results of this article and describe the global shapes of attainable 
yield curves in any affine one-factor term structure model. 

Definition 3.6. The (zero-coupon) yield Y(r t ,x) is given by Y(r t ,0) := r t and 

, qs x vf , \ogP{t,t + x) A{x) B(x) 

(3.8) Y(r t ,x) := = r t for all x > . 

x xx 

For r t fixed, we call the function Y(r t , .) the yield curve. 

The (instantaneous) forward rate f(r t , x) is given by f(r t , 0) := r t and 

(3.9) f(r t ,x) :=-d x log P(t, t + x) =- A' (x)-r t B'(x) for all x > . 
For rt fixed, we call the function f(rt, ■) the forward rate curve. 
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By l'Hospital's rule and the generalized Riccati equations ()3.6j) it is seen that 
both the yield and the forward rate curve are continuous at 0. 
The first quantity associated to the yield curve that we consider, is the asymptotic 
level frasymp of the yield curve as x — > oo, also known as long-term yield, consol 
yield or simply 'long end'. 

Theorem 3.7. Let the short rate process be given by a one- dimensional affine pro- 
cess (r t )t>a satisfying Condition \S.1\ with quasi-mean-reversion A. 
If X > then 

&asymp := lim Y(r t ,x) = lim f(r t ,x) = -F(-l/X) . 

X^OO X — >oo 

If X = then 

&asymp = hm -F(u) + T t (1 - R(u)) . 
u — > — oc 

Proof. From (|3.6ap we obtain that 

A(x) 

(3.10) lim — — = lim A'{x) = lim F(B(x)) . 

x — ► oo x x— »oo x — *oo 

If A > then by Corollary [33] 

B(x) 

(3.11) lim B{x) = -1/A, lim — >-L = and lkB'(i)=0 

x — >oc x — >oo X x — >oo 

and the assertion follows by combining l|3.8p - l|3.1ip . 
If A = then lim^^oo B(x) — — oo and 

li m EM. = ii m B'(x) = lim R(B{x)) - 1 . 

x — >oo x x — *°° x — 

By setting u := B(x) we obtain the desired result. □ 

From Theorem 13. 71 it is clear that for practical purposes only models with A > 
will be useful. So far we know that in this case the short end of the yield curve is 
given by Y{r t , 0) — r t and the long end by Y(r t , oo) = & a symp- We will now examine 
what happens between these two endpoints. 

Definition 3.8. The yield curve Y(r t ,x) is called 

• normal if it is a strictly increasing function of x, 

• inverse if it is a strictly decreasing function of x, 

• humped if it has exactly one local maximum and no minimum on (0,oo). 
In addition we call the yield curve flat if it is constant over all x € M^o- 

This is our main result on the shapes of yield curves in affine one-factor models: 

Theorem 3.9. Let the risk-neutral short rate process be given by a one- dimensional 
affine process (rt)t>o satisfying Condition ] 3. H and with quasi-mean-reversion A > 0. 
In addition suppose that F ^ and that either F or R is non-linear. Then the 
following holds: 

• The yield curve Y(r t , .) can only be normal, inverse or humped. 



• Define 



F' -1/A , , i-lM W° <° 

— — — - and b inv := < R'(0) 

K y- 1 ' X > [+O0 ifR'{0)>0 
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The yield curve is normal if rt < frnorm , humped if 6 n0 rm < r* < b mv and 
inverse if r t > b inv ■ 

The above theorem is visualized in Figure [TJ For its proof we will use a simple 
Lemma. We state the Lemma without proof, since it follows in an elementary way 
from the usual definition of a convex function on R. 

Lemma 3.10. A strictly convex or a strictly concave function on R intersects an 
affine function in at most two points. In the case of two intersection points p\ < pi, 
the convex function lies strictly below the affine function on the interval (pi,P2); if 
the function is concave it lies strictly above the affine function on (p\,p2). 

Proof of Theorem\3M Define the function H(x) : R^ -> R by 

(3.12) H(x) := Y{r t ,x)x = -A{x) - r t B(x) . 

We will see that the convexity behavior of H will be crucial for the shape of the 
yield curve Y(r tl .). From the generalized Riccati equations (|3.6p the first derivative 
of H is calculated as 

(3.13) d x H{x) = -F(B(x)) - r t (R(B(x)) - 1) 
and the second as 

(3.14) d xx H{x) = -B'{x) {F'{B{x)) + r t R'{B(x))) . 

Note that F and R are continuously differentiable by Lemma [2.6i and also B by 
(|3.6bj) . such that the second derivative of H is well-defined and continuous. Since B 
is strictly decreasing by Corollarv l3.5( the factor —B'(x) is positive for all x £ R^o- 
The sign of d xx H(x) therefore equals the sign of 

(3.15) k(x) := F'(B{x)) + r t R'{B{x)) . 

From the fact that B is decreasing and F and R are convex it is obvious that k 
must be decreasing. We will now show that k has at most a single zero in [0, oo): 

(a) D = K^o : We have assumed that either F or R is non-linear. By Lemma [231 
this implies that either F or R is strictly convex, and thus that cither F' or R' 
is strictly increasing. If rt > 0, then it follows that k is strictly decreasing and 
thus has at most a single zero. If r t = 0, an additional argument is needed: 
It could happen that F is of the form F = bu such that k(x) — b and k is no 
longer strictly decreasing. However, by assumption, F / such that in this 
case k has no zero in [0, oo). 

(b) D = R: In this case, by the admissibility conditions in Definition 12.31 we 
have necessarily R(u) — (3u. Also, since either F or R is non-linear, F must 
be non-linear and thus by Lemma 12.61 strictly convex. It follows that k(x) = 
F' (B(x))+rtf3 is strictly decreasing and thus has at most a single zero in [0, oo). 

We have shown that k is decreasing and has at most a single zero; to determine 
whether it has a zero for some value of rt, we consider the two 'endpoints' fc(0) and 
lim^^oo k{x). First we show that 

f ^(0) 

(3.16) fc(0) > if and only if r t < b inv := I R'(0) 

[+oo ifi?'(0)>0 

Since B(Q) — by Proposition 13.41 it follows that 

k(0) = F'(0) + r t R'(0) . 



if i?'(0) < 
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We distinguish two cases: 

(a) If i?'(0) < then the assertion (|3.16[) follows immediately. 

(b) Consider the case that R'(0) > 0: Assume that D = R. Then we have 
R(u) — j3u and R'(0) = f3 > 0. This, however, stands in contradiction to 
our assumption A > 0, which implies that (3 = — A < (cf. Definition 13. 2p . 
Thus we must have D = M.^ and r t > 0; in this case it follows that A;(0) > 0, 
for all r t € D, and we set 6 inv = +oo. 

Next we consider the right end of k(x) and show that 

F'(— 1/X) 

(3.17) lim k(x) < if and only if r t > b noim := - / . 

Since lim^^oo B(x) = — 1/A by Corollary 13. 51 we have that 

(3.18) lim k(x) = F'(-l/A) + r t R'(-l/X) . 

X — >oc 

By assumption A > 0, and by Definition 13.21 it holds that R{— 1/A) = 1. Also 
i?(0) = 0, and by the mean value theorem 

1 = R(-1/X)-R(Q) = ~R'(0 

for some £ e (— 1/A,0). Since R' is increasing, it follows that R'(— 1/A) < —A < 0, 
and we can deduce (|3.17j) directly from (|3.18| . 

We summarize our results on the function k so far: k stays negative on (0,oo) if 
rt > hnv and positive if rt < &norm- It has a single zero on (0, oo) if and only if 
&norm < r t < b mv . If k has a zero on (0, oo), since k is decreasing, the sign of k will 
be positive to the left of the zero and negative to the right of the zero. 

Since d xx H has the same sign as k, the statements above translate in the obvious 
way to the convexity behavior of H. We will now use the convexity behavior of H 
to derive our results about the yield curve. 
Consider the equation 

(3.19) H(x)=cx, xe[0,oo) 

for some fixed c € M. Since H(0) = this equation has at least one solution, 
.To = 0. If r t > &inv then H(x) is strictly concave on [0,oo), and by Lemma 13.101 
the equation (|3. 19|) has at most one additional solution x\. Also, when the solution 
exists, H(x) crosses cx from above at x\. Similarly if rt < frnorm then H(x) is 
strictly convex, and there exists at most one additional solution X2 to (|3.19p on 
[0, oo). If the solution exists, then cx is crossed from below at X2- In the last case 
frnorm < Tt < foinv, there exists a cc* - the zero of k(x) - such that H(x) is strictly 
convex on (0, x*) and strictly concave on (x*, oo). Now there can exist at most two 
additional solutions x±, x-i to (|3.19[) with x\ < x* < X2, such that cx is crossed from 
below at X\ and from above at X2- 

Because of definition (|3.12j) . every solution to (|3.19[) , excluding xq = 0, is also a 
solution to 

(3.20) Y(r t ,x) = c, x e (0,oo) 

with r t fixed. Also the properties of crossing from above/below are preserved since 
x is positive. This means that in the case r* > &inv, equation (|3.20p has at most a 
single solution, or in other words, that every horizontal line is crossed by the yield 
curve at most in a single point. If it is crossed, it is crossed from above. This 
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implies that Y(x) is a strictly decreasing function of x, or following Definition 13. 8[ 
that the yield curve is inverse. In the case rt < &norm we have again that (I3.20|l has 
at most a single solution and that every horizontal line is crossed from below by 
the yield curve, if it is crossed. In other words, the yield curve is normal. In the 
last case of b norm < r t < bi nv , the yield curve crosses every horizontal line at most 
twice, in which case it crosses first from below, then from above. Thus in this case 
the yield curve is humped. □ 



Corollary 3.11. Under the conditions of Theorem \3.9\ the instantaneous forward 
rate curve has the same global behavior as the yield curve, i.e. 

Y(r t , .) is inverse f(r t , .) is strictly decreasing 

. . , , f(rv,.) has exactly one local maximum 

Y(r t , .) is humped <^> , , , ■ 

and no local minimum 

Y(r t , .) is normal /(?'t, •) is strictly increasing . 

In the second case the maximum of the forward rate curve is f(rt,x*), where 
solves 

<3 - 21 > r ' = -wm- " (0 ' oo) - 

Proof. This follows from the fact that d x H{x) as given in (|3. 13|) is exactly the 
forward rate f{r t ,x). The derivative of the forward rate is therefore d xx H(x), 
which is given in (|3.14p as 

d x f(r t ,x) = d xx H(x) = -B'(x) ■ k(x) . 

The factor —B'(x) ^ is always positive, and the possible sign changes and zeroes 
of k(x) are discussed in the proof of Theorem l3.91 leading to the stated equivalences. 
Equation (|3.21|) is simply the condition k(x*) =0. □ 



Corollary 3.12. Under the conditions of Theorem \3.9\ it holds that 

(3.22) ^norm ^ ^asymp ^ ^inv 

whenever the quantities are finite. In addition it holds that 

(3.23) D ("I (&norm, Vv) £ • 

Remark 3.13. Note that equation (13.23|) implies that there is always some r t G D 
such that the yield curve Y(r t , .) is humped. 

Proof. By the mean value theorem there exists a £ S (— 1/A, 0) such that 

&a S ymp = ~F(-1/X) = F(0) - F(-l/A) = \f> (£) . 

Since F is convex and thus F' increasing, it holds that 

, nn * F'i-l/X) , F'(0) 

(3.24) \ ~ 6asymp " — \ ■ 

Applying the mean value theorem to i?, there exists another £ G (— 1/A, 0) such 
that 

l = R(-l/\)-R(0) = —#($). 
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Since R! is increasing we deduce that R'(—l/X) < —A < 0. Assuming that also 
R'(0) < we get 

(3 ' 25) " R'(-l/X) " A " ~R(0) ' 

Since either F or R is non-linear, one of the functions is strictly convex by Lemma [2.6l 
Consequently either both inequalities in (|3.24[) or in (|3.25p are strict. Putting them 
together we get 

F'(-l/X) F'(0) 
R'(-l/\) < ° asymp < R'(Q) ' 
proving (I3.22|) under the assumption that R'(0) < 0. 

If R'(0) > then by definition 6i nv = oo. Equation (I3.24| still holds, but in 
(|3.25p only the left inequality sign remains valid. Together this still proves that 
frnorm < &asymp and we have shown (|3.22[) . 

To prove (|3.23|) we distinguish two cases: 

(a) D = M. In this case it is sufficient to prove — oo < b- mv and b aotm < oo. Consider 
first &inv If R'{0) > then by definition b mv — oo and nothing is to prove. 
If R'(0) < then 6 inv = -F'(0)/R'(0). By convexity F'(0) > -oo and the 
assertion follows. Consider now b noim = —F'(—l/X)/R'(—l/X). From (|3.25|) 
we know that R'(— 1/A) < —A < 0. By convexity F'(—l/X) < oo and it follows 
that 6 norm < oo. 

(b) D = M^o- In this case it is sufficient to prove < 6 norm and to apply (|3.22|) . 
As above we have that 6 norm = -F'(-l/X)/R'(-l/X) and that R'(-l/X) < 
-A < 0. By Definition EH 

F'(-1/X) = b+ I £e-& x m(d4) 

with b > 0. It follows that F'(—l/X) > 0, proving the assertion. □ 

The last Corollary of this section shows the interesting fact that the occurrence 
of a humped yield curve is a necessary and sufficient sign of randomness in the 
short rate model: 

Corollary 3.14. Let the risk-neutral short rate process be given by a one- dimensional 
affine process (rt)t>o satisfying Condition \3.1\ with F ^ and quasi-mean-reversion 
A > 0. Then the following statements are equivalent: 

(i) There exists a r t £ D such that Y(r t , .) is flat. 

(ii) There exists no r t S D such that Y(r tl .) is humped. 
(Hi) The short rate process (r t )t>o is deterministic. 

(iv) F{u) = bu and R(u) — (3u. 

Proof. Theorem l3.91 together with Corollarv l3.121 shows already that -i(iu) implies 
->(i) and Also, from the form of the generator in (|2.4p . it is seen that (Hi) and 

(iv) are equivalent. It remains to show that (iv) implies (i) and (ii). Proceeding 
as in the proof of Theorem 13.91 we obtain instead of (|3.15[) simply 

k(x) = b + r t (3 . 

The yield curve will be humped if and only if k has a single (isolated) zero in [0, oo). 
Since A: is a constant function, this cannot be the case for any r* € D and we have 
shown (ii). By the same arguments as in the proof of Theorem 13.91 the yield curve 
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FIGURE 1 . This Figure shows a graphical summary of Theorems 
13.71 and 13. 9( as well as the definitions of the key quantities b norm , 
feaaymp and 6j nv . In any affine model satisfying the conditions of 
Theorem 13.91 the shapes of yield curves will follow the picture 
given here. They will be normal if ro is below 6 norm , humped if r$ 
is between 6 n orm and b- lnv and inverse if ro is above b lnv . Also all 
yield curves will tend asymptotically to the same level & a symp- 



is fiat if and only if k is constant and equal to 0. This is the case if 7"t = — 4. It 
remains to show that rt € D. Note that (3 = — A < 0. In particular /3 ^ 0, such 
that for D = R we are already done. If D = R^o we have by the admissibility 
conditions in Definition 12.31 that b > 0. Thus rj = — \ > and we have shown 
(i). ' □ 



3.3. The Limit Distribution of an Affine Process. It is well-known that the 
Gaussian Ornstein-Uhlenbeck process, for example, converges in law to a limit 
distribution and that this distribution is Gaussian. The goal of this section is to 
establish a corresponding result for affine processes. While calculating the marginal 
distributions of an affine process involves solving the generalized Riccati equations 
(|2.5[) . it will be seen that the limit distribution is much easier obtained and can be 
determined directly from the functions F and R. 

In the interest rate model considered in the preceding section, the short rate follows 
an affine process under the martingale measure, such that the results will allow us 
to characterize the risk-neutral asymptotic short rate distribution. Often also the 
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limit distribution under the objective measure is of interest, but the affine prop- 
erty is in general not preserved by an equivalent change of measure, such that the 
results are not directly applicable. Nevertheless, for the sake of tractability, condi- 
tions on the measure change can be imposed, such t hat the model is affine under 
both the objective and the risk-neutral measure. (SeelNicolatq and Venardosl [2003j 



for an example from option pricing and ICheridito et al. 20051 ] for more general re- 



sults). In such a setting the results can also be applied under the objective measure. 

Before we state the result, we want to recall that a real- valued random variable 
L is called self-decomposable if for every c £ (0,1) there exists a random variable 
L c , independent of L, such that 

L = cL + L c for all c e (0, 1) . 

Since self-decomposability is a distributional property, we will identify L and its 
law, and refer to both as self-decomposable. 

For OU-type processes, limi t distributions have been stu died for some time; 
the fir st results can be found in Jurek and Vervaat 1983j and Sato and Yamazatol 



1984]. The next t heorem summarizes these results, and can be found in similar 



form in lSatd [1999L Theorem 17.5] 



Theorem 3.15. Let (rt)t>o be a OU-type process on R. // 



P < and / log |£| m{dg) < oo 
J\S\>1 

then (rt)t>o converges in law to a limit distribution L which is independent of ro 
and has the following properties: 
(i) L is self- decomposable. 

(ii) The cumulant generating function k(u) — log J R e ux dL(x) satisfies 

1 f° F(is) 

(3.26) K(iu) = - / — —ds for all u£l. 

P Ju s 

Conversely, if L is a self-decomposable distribution on R and [3 < 0, there exists 
a unique triplet (a, 6, m) satisfying the admissibility conditions of Definition \2.3l 
such that L is the limit distribution of the affine process (of OU-type) given by the 
parameters (a, 6, m, /?). 

As discussed in Section [21 every regular affine process with state space R is of 
OU-type, such that the above theorem applies. We now state our corresponding 
result for affine processes on R^o: 

Theorem 3.16. Let (rt)t>o be a one- dimensional, regular, conservative affine pro- 
cess with state space R^o- If 



R'(0) < and / log£m(d£) < oo 
J£>1 



then (rt)t>o converges in law to a limit distribution L which is independent of ro, 
and whose cumulant generating function n is given by 

f° F(s) 

(3.27) k(u)= -44 ds for all ug(-oo,0]. 

Ju 
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Proof. By Theorem 12.41 the transition kernel pt(x, .) of the process (n)t>o has the 
characteristic function 

Pt{x, u) = exp (4>(t, u) + xt{)(t, u)) 

where </> and ip satisfy the generalized Riccati equations (|2.5j) for all u E U, and 
thus in particular for all u E (— oo,0]. Since R(0) = 0, is a critical point of the 
autonomous ODE (I2.5bj) . and by the assumption i?'(0) < it is asymptotically 
stable. By the convexity of R, R'(Q) < also implies that R(u) > for all u E 
(-co, 0), such that tp(t, u) is a strictly increasing function in t for all u E (— oo, 0). 
Since is the only critical point of (|2.5bl) on (— oo,0] it also follows that 

lim ip(t, u) = for all u E (— oo, 0] . 

t — >oc 

Consequently, 

f°° f° F(s) 

(3.28) lim logp t (x,u) = lim (f>(t,u) = / F(ip(r,u))dr = —^4 ds 

t-oo t-oo J Q J u R(s) 

where the last two equalities follow from l|2.5p and the transformation s = ip(r,u). 
We will now show that the last integral in (|3.28p converges absolutely for all u E 
(— oo, 0]: Since R(u) > and F(u) < for all u E (— oo, 0] we obtain 

° F ( s ) ds <__L_ /°£Md Si «G (_oo,oi, 



/; 


F(s) 


ds = - [ 


/ u 


R(s) 


J U 



R(s) ~ R'(0) 

where the inequality follows from the fact that the convex function R is supported by 
its tangent at 0. From the definition of F(u) in (|2.2p it is clear that the convergence 
of the last integral will depend only on the jump part of F, i.e. the integral converges 
if and only if 

(3.29) /-/ ( e ' Si ~ l)m{d£)ds < oo, for all u G (-oo, 0]. 

in s 7(0, oo) 

Define M[u, £) = f° e ( ~ 1 ds. For a fixed u G (— oo,0], it is easily verified that 
M(u,g) = 0(0 as £-> 0, and that M(u,£) = O(log^) as £ -> oo. Since the Levy 
measure m(d£) integrates (£ A 1) by Definition ^. 3i and log^- 1{£>i} by assumption, 
it must also integrate M(u, £). Applying Fubini's theorem, (|3.29[) follows, such that 
k(u) := J u S|t ds converges for all u E (— oo, 0]. In particular lim„fo k(u) — 0, such 
that the limit in (|3.28|) is a function that is left-contin uous at 0. By standard r esults 
on Laplace transforms of probability measures (cf. ISteutel and van Harn 2004 . 



Theorem A. 3.1]), the pointwise convergence of cumulant generating functions to a 
function that is left-continuous at implies convergence in distribution of (r t ) t >o 
to a limit distribution L with cumulant generating function given by (13.28|) . □ 

Since the marginal distributions of an affine process are infinitely divisible, also 
the limit distribution L must be infinitely divisible, if it exists. In Theorem 13.151 
a stronger result is given for an affine process on M: In this case L is also self- 
decomposable. An obvious question is, if this result can be extended to the state 
space Mj>o- We will see that the answer is negative. In Section [4~3l an example of an 
affine process with state space IR^o is given, which converges to an infinitely divisible 
limit distribution that is not self-decomposable. This result is interesting, since it 
leaves open the possibility of some unexpected properties of the limit distribution of 
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an affine process. For example a self-decomposable distribution is always unimodal, 
whereas an infinitely divisible distribution might be not. 

4. Applications 

4.1. The Vasicek model. We apply the results of the last section to the classical 
Vasicek model 

(4.1) dr t = -\(r t -6)dt + adW t , r e K 

where (Wt)t>o is a standard Brownian motion under the risk- neutral measure and 
A, 9, a > 0. The Vasicek model is arguably the simplest affine model, and no 
surprises are to be expected here . In fact all re sults that we state here can already 
be found in the original paper of Vasicek! 1977 1. We advise the reader to view this 



paragraph as a warm-up for the following examples. 
Clearly (r t ) t >o is a conservative affine process with 

2 

(4.2) F(u) = X9u + y« 2 , 

(4.3) R(u) = -Au . 

From the quadratic term in F and Definition 12.31 it is seen that (rt)t>o has state 
space HL This property is often criticized, since it allows the short rate to become 
negative. 

From Theorem 13. 9l we calculate 

2 

hnv = @ and fonorm = — -r~ , 

such that the yield curve in the Vasicek model is normal if r t < 9 — er 2 /A 2 , inverse 
if Tt > 9 and humped in the remaining cases. 
The long term yield is calculated from f|3 . T|) as 

^_ 

2A 2 ' 

in this case exactly the arithmetic mean of 6; nv and ii norm . 

Theorem 13.151 applies and the cumulant generating function k of the risk-neutral 
limit distribution L satisfies 

^4/;^-i>-£»)--- % 

for Hence, L is Gaussian with mean 6 and variance Sy- 



&asymp — — F{— 1/A) 



4.2. The Cox-Ingersoll-Ross model. The Cox-Ingersoll-Ross (CIR)-model was 
introduced by lCox et al.l [1985J . In this model the short rate process (rt)t>o is given 
by the SDE 

(4.4) dr t = ~a{r t - 6)dt + a^r~ t dW t , r € M >0 

where (Wt)t>o is a standard Brownian Motion under the risk-neutral measure and 
a, 9, a > 0. The process (r t ) t >a is a conservative affine process with 

(4.5) F(u) = a9u , 

2 

(4.6) R{u) = - au . 
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From Definition 12.31 it is seen that {rt)t>o has state space K^o- The fact that 
interest rates stay non-negative in the CIR-model is often cited as an advantage 
of the model over the Vasicek model. Calculating the quasi-mean-reversion (see 
Definition 13. 2|) . we find that 



From Theorem 13 . 71 we find that the long-term yield is given by 

2a9 

Oasymp — — F(— 1/A) — 



\Ja? + 2cr 2 + a ' 

The boundary between humped and inverse behavior b lnv is calculated from Theo- 
rem 13.91 as 

inv " R'(0) - 

Both quantities b asymp and b- mv can also be found in [Cox et al. , 1985. Eq. (2G) and 



following paragraph]. Befo re we cons i der 6n ormi we quote (with notation adapted 
to (|4.4[0 from page 394 of Cox et all Il985j where the shape of the yield curve is 



discussed: 

'When the spot rate is below the long-term yield [= 6 asymp ], the 
term structure is uniformly rising. With an interest rate in excess 
of 9 [= &in V ], the term structure is falling. For intermediate values 
of the interest rate, the yield curve is humped. ' 

In our terminology, they claim that the yield curve is normal for r t < & asymp , 
humped for & asymp < r t < 6; nv and inverse for r t > b inv . This stands in clear 
contradiction to Theorem 13 .91 and Corollarv l3 . 1 21 where we have obtained that yield 
curves are normal if and only if r t < b noim and that b noim < fe asym p, or - in plain 
words - that there are yield curves starting strictly below the long-term yield that 
are still hum ped. 

The claims o f lCox et all [1985( | are repeated in |RebonatoL Il998l p. 244f], where even 



several plots of 'yie ld surfaces ' (the yield as a function of r t and x) are presented as 
evidence. However iRebonatol fails to indicate the level of & asym p in the plots, such 
that the conclusion remains ambiguous. 

To clarify the scope of humped yield curves in the CIR-model we calculate 6 n orm 
from Theorem 13.91 

F'j-l/X) _ a9 

Onorm- # ( _ 1/A) ~ ' 

The relation 6 norm < & asym p < ^inv is immediately confirmed by noting that 6 asymp 
is the harmonic mean of 6 norm and b- mv . For a graphical illustration we refer to 
the second yield curve from below in Figure [1] The plot actually shows CIR yield 
curves with parameters 

a = 0.5, cr = 0.5, 9 = 6% 



plotted over a time scale of 25 years. The second curve from below starts at 
tq = 4.2%, i.e. below the long-term yield, but is visibly humped. 
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To calculate the limit distribution of (rt)t>o, we apply Theorem 13. 161 The cu- 
mulant generating function k(u) of the limit distribution is given by 

f° F ( s ) , r , 2ad ( a 2 
< u ) = / "57^ ds = / i ds = ZT lo § 1 ~ 7T7 U 



R(s) J l-sa 2 /2a a 2 b \ 2a 

This is the cumulant generating function of a gamma distribution with shape pa- 
ra meter 2a9/a 2 an d scale parameter a 2 /2a. Again this result can already be found 
in lCoxet alJ |l985l p. 392]. 



4.3. An extension of the CIR model. To illustrate the power of the affine 
setting, we consider now an extension of the CIR model that is obtained by adding 
jumps to (|4.4p . We define the risk-neutral short rate process by 



(4.7) dr t = -a(r t - 6)dt + a^/n dW t + dJ t , r > 

where (Jt)t>o is a compound Poisson process with intensity c > and expo- 
nentially distributed jump s of mean v > 0. This model has been intro duced by 



Duffle and Garleanul [200 lj as a model for default intens ity and is used by Filipovic 
200l| as a short rate model. It can also be found in Bri go and Mercurid [2006 



under the name JCIR model. It is easily calculated that 

cu 

(4.8) F(u) = a0u-\ , u6(-oo,v), 

v — u 

2 

(4.9) R(u) = yii 2 - au . 

Solving the generalized Riccati equations (|3.6[) for A{x) and B(x) becomes quite 
tedious, but the quantities b- mv , 6 a symp , b norm can be calculated from Theorem 13.71 
and Theorem 13.91 in a few lines: The quasi-mean reversion A stays the same as in 
the CIR model, since R does not change. From 



F'(u) = a6 + 



{v-u) 2 



we derive immediately 



b, 



binv — 9 H , 

av 

2a0 2c 



bnnrm — 



a + 7 v{a + v) + 2 ' 
a9 



4 

cva 



7 j(a 2 v + 7 — a) 2 ' 

where 7 = \J a 2 + 2a 2 . Note that by setting the jump intensity c to zero, the ex- 
pressions of the (original) CIR model are recovered. 

Next we calculate the limit distribution of the model. Using the abbreviations 
p := a 2 /2 and A := a — vp we obtain 

« U) = [°^d S = r°d S .-~r ds 



u R{s) Jo l-sp/a J {s-v)(ps-a) 



A p 



log (1 log if A^0 

-0„log(l-a) + £_JL. if A = 
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as the cumulant generating function of the limit distribution L under the martin- 
gale measure. 

We now take a closer look at the distribution L, since this will answer the 
question raised at the end of Section [3.31 For certain parameters, L is an example 
for a limit distribution of an affine process that is infinitely divisible, but not self- 
decomposable. We consider the case A = and define 

(4.10) l{x) := (6 + -x) ve' vx , x G 



a 



By Frullani's integral formula 

r°° l(x) 
(4.11) k(u)= / (e vx -i) l -¥±dx 

' x 







for all u G (—00,1/). Since I is non-negative on K^cb l{x)/x is the density of a 
Levy measure and (|4.1ip is seen to be the Levy-Khintchine representation for the 
cumulant generating function of the infinitely divisible distribution L. In addition, 
L is self-decom posable if and only if I is non-negative and non-increasing on R^ 
(cf. ISatol (1999L Corollary 15.11]). 



In the case of l(x) given by (|4.10[) . it is easily calculated that l(x) has a single 
maximum at i* = - — — . Thus, if c < aOv, then x* < 0, such that I is non- 

v c ' — 5 — 5 

increasing on R^o an d L is self-decomposable. If c > a8v then I is increasing 
in the interval [0, x*) and the limit distribution L is infinitely divisible, but not 
self-decomposable. 

4.4. The gamma model. Instead of analyzing the properties of a known model, 
we will now follow a different route and construct a model that satisfies some given 
properties. We want to construct an affine process on R^o that has the same limit 
distribution as the CIR model (i.e. a gamma distribution), but is a process of OU- 
type. The second property is equivalent to R(u) = [3u. Considering Thcorcm l3.161 
we know that if we want to obtain a limit distribution, we need < 0. To keep 
with the notation of the Vasicek model, we will write R(u) = — Xu where A > 0. 
Now by (|3.27j) the cumulant generating function of the limit distribution is given 

by 

1 f u F(s) 

(4.12) k(u) = - / — —ds for all u G (-oo,0l . 

A Jo s 

Let the limit distribution be a gamma distribution with shape parameter k > 
and scale parameter 8 > 0. Then k(u) — — fclog(l — 6u) and by (|4.12p 

\9ku 

F{u) = T3fe ' 

Setting c = Xk and v = 1/8 it is seen that F(u) is equal to the last term in 
(|4.8[) . This means that the driving Levy process of (r t ) t >o is of the same kind 
as the process (Jt)t>o in (|4.7|) . i.e. (r t )t>o is a pure jump OU-type process with 
exponentially distributed jump heights of mean 1/8 and with jump intensity Xk. 
We interpret the affine process we have constructed as a risk-neutral short rate 
process. It is clear that the bond prices are of the exponentially-affine form 
From the generalized Riccati equation (|3.6b[) we obtain 



B(x) 



e -\x _ l 

A 
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From equation (|3.6a|) wc calculate 

f x \k 
A{x) = J F(B(s)) ds = — (log(l - 0B{x)) - 9x) , 

such that the bond prices are given by 

P(t,t + x) ^exp!^-x-^j+r t B(x)^ {l-eB(x))T^ . 

The global shape of the yield curve is described by the quantities 

k k/0 

6i„v - k9, 6 asymp - 1/9+1/ y knonn - QfQ+yXf 

and it is seen that for the gamma-OU-process fo a symp is the geometric average of 

^inv and fenorm ■ 



5. Conclusions 

In this article we have given, under very general conditions, a characterization of 
the yield curve shapes that are attainable in term structure models where the risk- 
neutral short rate is given by a time-homogenous, one-dimensional affine process. 
Even though the parameter space for this class of models is infinite-dimensional, 
the scope of attainable yield curves is very narrow, with only three possible global 
shapes. In addition we have given conditions under which an affine process con- 
verges to a limit distribution, and we have characterized the limit distribution in 
terms of its cumulant generating function, extending some known results on OU- 
type processes. 

The most obvious question for future research is the extension of these results to 
multi- factor models. It is evident from numerical results that in two- factor models 
yield curves with e.g. a dip, or also with a dip and a hump, can be obtained. It 
would be interesting to see if more complex shapes can also be produced, or if there 
are similar limitations as in the single-factor case. Also, in the one-factor case the 
dependence of the yield curve shape on the current short rate is basically described 
by the intervals D fl (— oo, 6 n orm], (&norm, ^inv) and [b- mv , oo). In the two-factor case 
the partitioning of the state-space might be more complex, and we expect to see 
more interesting transitions between yield curve types. Another aspect is, that 
since affine processes as a general framework become better understood, extensions 
of classical models e.g. by adding jumps, like in the JCIR model described in 
Section [4.31 become more feasible and attractive for applications. 
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